Nonlinear stochastic discrete drift-diffusion theory of charge fluctuations and domain 

relocation times in semiconductor superlattices 
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A stochastic discrete drift-diffusion model is proposed to account for the effects of shot noise 
in weakly coupled, highly doped semiconductor superlattices. Their current-voltage characteristics 
consist of a number stable multistable branches corresponding to electric field profiles displaying 
two domains separated by a domain wall. If the initial state corresponds to a voltage on the middle 
of a stable branch and a sudden voltage is switched so that the final voltage corresponds to the next 
branch, the domains relocate after a certain delay time. Shot noise causes the distribution of delay 
times to change from a Gaussian to a first passage time distribution as the final voltage approaches 
that of the end of the first current branch. These results agree qualitatively with experiments by 
Rogozia et al (Phys. Rev. B 64, 041308(R) (2001) ). 
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I. INTRODUCTION 



The current-voltage (I-V) characteristics of highly 
doped weakly coupled semiconductor superlattices (SL) 
typically exhibits many sharp branches due to formation 
of static electric field domains ^ . Two branches are sep- 
arated by a discontinuity in the current. The electric 
field profile associated to a given branch consists of two 
regions of constant electric field (domains) separated by 
a charge accumulation layer (domain boundary), which 
is confined to one or several quantum wells. The location 
of the domain boundary distinguishes I-V branches: as 
the voltage increases, the domain boundary is located 
closer to the injecting contact and the high field do- 
main increases at the expense of the low field one 
Branches exhibit hysteresis cycles due to coexistence of 
two or more stable electric field profiles at a given value of 
the voltage. Many interesting dynamical phenomena are 
associated to these SL: (i) response of the SL to sudden 
changes in bias (which may force relocation of electric 
field domains |^-^ ) , and (ii) self-sustained oscillations of 
the current provided temperature is raised or doping is 
lowered [|^,^ . Motivated by recent experimental evidence 
[^,|o|, we shall present in this paper a stochastic theory 
of domain relocation in highly doped SL. 

In relocation experiments a doped SL dis- 

playing a multistable I-V characteristic is biased (typ- 
ically) on the first plateau, say in the middle of a branch. 
The corresponding field configuration has two domains 
separated by a domain wall which is an accumulation 
layer. Then the voltage is suddently increased from Vq 
to Vi = Vo-\- AV and the time evolution of the current is 
recorded. Depending on AT^, the domain wall has to relo- 
cate so that a stable field configuration appropriate to the 



new voltage is reached The outcome has been stud- 
ied numerically using a discrete resonant tunneling model 
with Ohmic boundary conditions For any AV < 
as well as for small positive AV, the relocation of the 
domain wall always occurs by a direct movement of the 
charge monopole forming the domain boundary to its fi- 
nal position. This movement may be either upstream or 
downstream the electron flow as needed. However, for 
sufficiently large AV > 0, a charge dipole is injected at 
the emitter contact in addition to the existing monopole, 
because the latter cannot move upstream beyond one SL 
period without encountering a stable field configuration 
Recent experiments by Rogozia et al. jl^] confirm 
this theoretical picture. Other experiments have shown 
that the relocation time for up jumps {AV > 0) close to 
the discontinuity in the I-V characteristic is random and 
have also investigated its probability distribution func- 
tion |9[[Tc[| . What is causing randomness in the relocation 
time? In this paper we argue in favor of shot noise. 

Shot noise occurring during a transport process is 
due to fluctuations in the occupation number of states 
caused by (i) thermal random initial fluctuations; (ii) 
the random nature of quantum-mechanical transmis- 
sion/reflection (partition noise). The latter is in turn 
caused by the discrete nature of the electric charge. 

The rest of the paper is organized as follows. In Section 
H, we derive a stochastic discrete drift-diffusion model 
(DDD) from the previously studied deterministic one (see 
Ref. |12| ) considering partition only noise (thermal noise 
is negligible in the low temperature limit) . The stochastic 
DDD model has multiplicative white noise terms obeying 
Poissonian statistics and it has been solved numerically 
by means of a second order scheme proposed by Platen 
p6[ . The results of numerically solving the stochastic 
model are reported in Section [II. Our numerical re- 
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suits agree qualitatively with Rogozia et al experiments 
1^ , thereby enforcing the idea that shot noise is respon- 
sible for the observed fluctuations in domain relocation 
time. Details on the numerical scheme and comparison 
to rougher schemes and to the results of solving the de- 
terministic model with random initial conditions arc con- 
tained in the Appendix. 



II. STOCHASTIC DISCRETE DRIFT-DIFFUSION 
MODEL 



In weakly coupled SL, typically the scattering times are 
much shorter than the escape times from quantum wells. 
In their turn, the latter are shorter than typical dielectric 
relaxation times. This implies that the dominant mech- 
anism of vertical charge transport is sequential resonant 
tunneling and that the tunneling current across barriers 
can be considered to be stationary. An appropriate dis- 
crete model consists of the Poisson and charge continuity 
equations for the two-dimensional electron density rii and 
average electric field Fi at each SL period M: 



= -(n,-iV^) i = l,---,7V, (1) 

e 

dTi 

= Ji-i^i - Ji^i+i , i = l,---,N. (2) 



Here N]^, e and eJ^^i+i are the 2D doping density at 
the ith well, the average permittivity of the SL and the 
tunneling current density across the ith barrier, respec- 
tively We can differentiate Eq. (|l|) with respect to 
time and eliminate rii by using Eq. (H). The result can 
be written as a form of Ampere's law for the balance of 
current 



e dFi 
e dt 
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FIG. 1. Drift velocity, diffusion coefficient. 

The DDD model given by Eqs. (|), (H) and .h^i+i = 

Ji^i+i has a conceptual difficulty coming from charge 
quantization that motivates the introduction of shot 
noise terms. The electric charge in each SL period, eriiA 
[A is the SL cross section), should be a multiple of the 
electron charge e. This implies that the true charge fluc- 
tuates about the mean value given by the deterministic 
DDD model. To analyze charge fluctuations, we may use 
the Langevin ideas and add an appropriate stochastic 



term to J], 



id) 



The SL cross section A is very large (a 



circular cross section of diameter 120 /im wide as com- 
pared to a SL period of Z = 13 nm) and the barrier trans- 
mission coefficient is very small. Then we may use the 
classic Poissonian shot noise to model charge fluctuations 

Ji^i+l - j h Ji^^+i[t), (4j 

(r) 

for i ~ 1, ■ ■ ■ N ~ 1, where Jj_,j_|_]^ represents the random 
current which satisfies 



+ Ji- 



Jit). 



(3) 



(4:2.+i> = 0, 



(5) 



Here eJ{t) is the total current density through the SL, 
equal for all SL periods, and edFi/dt is the displace- 
ment current at the ith SL period. To have a closed 
system of equations, we need a constitutive relation link- 
ing eJi—,i+i to the unknowns rii and Fi. For a weakly 
coupled SL the stationary sequential tunneling current 
has been calculated by the Transfer Matrix Hamiltonian 
method or by the Green function formalism [Q. 
In both cases, for sufficiently high temperature J^^i+i 
may be approximated by a discrete drift-diffusion law, 
= n,v{F,)/l - D{F,) {n,+i - n,)/P The SL 
period is I and the electron drift velocity v{F) and the 
diffusion coefficient D{Fi) have the forms depicted in Fig. 

0- 



6,,S{t - t') {Al)-^[n,v(f) (F,) + (F,)], (6) 

and v'-''\ v^^'^ are defined as follows 

yib)(^F) = ^^^, v^f'>{F)^v^'''>{F)+v{F), (7) 

The logic behind this form of the random tunneling cur- 
rent is as follows. We consider that uncorrelated elec- 
trons are arriving at the ith barrier with a distribu- 
tion function of time intervals between arrival times that 
is Poissonian Then the shot noise spectrum for 

the current eJi^i_^iA is given by the average current. 
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[nit;(-^)(Fj) + ni+ii'(*)(i^i)] e^AJl, which in turn yields Eq. 
@. As remarked in Ref. ||l5|, this procedure assumes 
low transmission through barriers and it yields an upper 
bound for the shot noise amplitude. In addition, the 
tunneling current is approximated by a discrete drift- 
diffusion expression whose transport coefficients (drift 
velocity, diffusivity, . . . ) will be quantitatively differ- 
ent from those of the actual sample used in experiments. 
Given the exponential dependence of several quantities, 
relatively small differences in the location of extrema of 
the drift velocity, etc. may produce substantial differ- 
ences. Thus, the mathematical model provides quanti- 
tative differences in the results but it yields the correct 
qualitative behavior. 

The special nature of the emitter and collector layers is 
considered in the boundary conditions, given by Eq. (|^) 
with i = and i ~ N and different constitutive relations 
for the tunneling currents p^: 



JV 
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, TV) 



J. 



N^N+l — 



I 



(8) 
(9) 



Here we still have (Ji'Lii^i) = for i = and i = N, 
while the correlations are: 



Al 



6{t-t'), (10) 



The emitter current density eje , the emitter backward 
velocity lu^^' and the collector forward velocity w'-'') are 
functions of the electric field depicted in Fig. || . 
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FIG. 2. Current-field contact characteristics. 



In addition to the boundary conditions, the Ampere 
and Poisson equations should be supplemented with the 
voltage bias condition. 



V. 



(12) 



where V denotes voltage. Eqs. (|l|), (|), (|) to (y|) form 
a closed system of stochastic equations for n^, Fi and 
J. They constitute the stochastic DDD model. To ana- 
lyze this model, it is convenient to render all equations 
dimensionless. Let us denote by {Fm,vm) the coor- 
dinates of the first positive maximum of the drift ve- 
locity v{F). We adopt Fm, N^, vm, vmI, gN^^vm/I 
and eFmI/ {gN'^vm) as units of Fi, Ui, v{F), D{F), eJ 
and t, respectively. In our numerical calculations, we 
have used parameters corresponding to the 9/4 SL of 
Ref. at a temperature of 5 K. The three first sub- 
bands have energies of 44, 180 and 410 meV, respec- 
tively, and we assume that the spectral functions of the 
wells are Lorentzians with half-widths of 10 mcV p2[ . 
Then we find Fm = 11.60 kV/cm, = 1.5 x 
cm~^, vm — 718 cm/s, vmI — 9.33 x lO^'^ cm^/s and 
cN'^vm/I = 13.27 A/cm^. For a circular sample with a 
diameter of 120 /xm, the units of current and time are 
1.501 mA and 1.021 ns, respectively Then Eqs. (0), (||), 
(0) to (111) become 



(13) 
(14) 
(15) 



Jit) 
J{t) 
Jit) 



dt 

+ay^niivt + Di) + DiUi+i ^i(i) 
dEo 
dt 



+ J,iE^)-W,iEQ)n^ 



+ay/JeiEo) + WeiEo) Ui ^oit) , 

dEN 



dt 



W,iEN)nN 



+a^yWciEN)nN ^Nit) 



1 ^ 



(16) 
(17) 



Here we have used the same symbols for dimensional and 
dimensionless quantities except for the electric field and 
the coefficient functions in the boundary conditions. The 
parameters u = eN^/ieFM) ~ 1.772, = V/iFuNl) 
and a = y^e/ieFMA) « 3.232 x lO"'^ are the dimen- 
sionless doping, the average electric field (bias) and the 
noise amplitude respectively, ^i(t) is a zero-mean Gaus- 
sian white noise with correlation (5i(i)^j it')) = 6ijSit—t') 

i^iit) = ^iitm)/V^, where the ^i(tm) are independent 
identically distributed (i.i.d.) normalized Gaussian ran- 
dom variables for each discrete time tm and At is the 
dimensionless time step). The rest of the coefficients in 
Eqs. (|l|) to (|l|) are defined by 



Vi = viEi) 



viFuE,) 



VM 
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JeiEo) 



DjFMEi) 

VmI ' 
ji^\FMEo)l 

N^VM ' 

W'^'-^FmEo) 
vm 

W'^f\FME N) 



(18) 



The previous system of equations can be further simpU- 
fied since the electron densities rii and the total current 
density J{t) can be expressed in terms of the electric 
field and the bias. Differentiating Eq. (^4|) with respect 
to time, and using Eqs. ([l^) and (16), we obtain an ex- 
pression for the total current density J{t): 

d(j) 1 dEi 1 . , ^ , ^ , 

dt N ^ dt TV ^ ^ + ^ 



nNWcjEN) 
N 



N-l 



^ \/niVi + [n^ + ni+i)Di^i{t) 



N 



y/nNWciEN)^Nit) . 



Then the total current can be written as 
J = Ji + J2 • 

4=1 



Jl 



dt 



N 



(J2)0 = 0, 

(■12)1 - 



N 



UNWcjEN) 

N 



l<i<N, 



(J2)W = 



ay/nNWc{EN) 
N 



(19) 



(20) 
(21) 

(22) 

(23) 
(24) 



We can now insert these equations in the Ampere equa- 
tions ( p^ to ( p^ ) and eliminate rii by using Eq. ( p^ ) 
thereby obtaining a stochastic differential equation of the 
following form: 



(25) 



for the {N + l)-dimensional vector electric field E = 
{Eo,...,EnY . Here S'(E) isa{N + l)x{N+ 1) matrix 
and H is a {N -I- l)-dimensional vector having obvious 
forms which we do not write explicitly for the sake of 
conciseness. 

The stochastic differential equation ( p5| ) has been nu- 
merically solved by using two different methods: a first 
order Heun scheme (modified Euler scheme) and the sec- 
ond order scheme proposed by Platen [n6|. The second 



numerical scheme is rather more costly, but we had to 
use it to avoid that numerical errors mask the effects 
due to charge fluctuations. Technical details on numeri- 
cal schemes and a comparison of their performances are 
given in Appendix The results of our simulations are 
reported in the next Section. 



III. NUMERICAL RESULTS 



We have numerically investigated the sample of Ref. 
that was used in the relocation experiments |^,|| . It con- 
sists of a TV = 40-period SL with 9-nm wide GaAs wells 
and 4-nm wide AlAs barriers, and 2D doping 7V]^ = 
1.5 X 10"^^ cm~^, at a temperature T = 5 K. We have 
solved numerically the nondimensional equations in the 
units and dimensionless parameters introduced in Section 

0- 
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FIG. 3. Static electric field profile at V=2.1 V. 



Figures § and ^ show a typical static electric field pro- 
file (with two coexisting domains) and the first plateau of 
the time averaged I-V characteristics (obtained by volt- 
age up sweeping). To ascertain the influence of charge 
fluctuations in domain relocation, we start by setting a 
stationary fleld configuration corresponding to a voltage 
Vb = 0.65 V on the lower branch of Fig. |4[ At time t = 0, 
the voltage increases (in one time step) to its final value 
Vf on the next I-V branch. 
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FIG. 4. Part of the first plateau of the I-V characteristics. 
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FIG. 5. Time trace of the current when the voltage is 
switched from Vo = 0.65 V to V> = 0.737 V. 



Time traces of the current are depicted in Figures |5| 
to |7[ Notice that the vertical scale has been augmented 
sufficiently to see the fluctuations of the current, that are 
typically about 0.02 in size. To compare our numerical 
results to experimental ones, we need to characterize the 
domain relocation times and their distribution function. 
After a voltage switch, each realization of the random so- 
lution of Eq. ( p5|) gives rise to jumps in the mean current 
as depicted in Figures I tog. We compare the time trace 
of the current (time averaged over intervals of five time 
dimensionless units) to the value of the current in static 
I-V branches. The first time to that the current time 
trace differs less than 5 x 10"'' dimensionless units from 
its final stationary value, we consider that the domain 
relocation has ended. The distribution of time delays 
to taken over many realizations is then recorded. For a 
large voltage switch, the time delay before the current 
falls from its initial value to its final level is shorter than 
for a smaller voltage switch; compare Figs. I and I The 
differences between the time delays involved in these two 
cases (about 40 ns) are smaller than those recorded in 
experiments |^ . These differences occur because of over- 
estimation of the field Fm and the shot noise amplitude 
by our theoretical calculations with respect to those of 
the experimental sample, as we mentioned before. 
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FIG. 6. As in Fig. U but with a final voltage Vf = 0.75 V. 
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FIG. 7. As in Fig. |, but with a final voltage Vf = 0.755 V. 



In Ref. 1^ it was claimed that the time delay depends 
exponentially on the difference between the final value of 
the stabilized current, /, and the maximum value of the 
current (or mimimum value in the case of a down switch) 
at the initial branch, /,„ . Then the relocation time (mea- 
sured in units of 1.021 ns) depends exponentially on the 
current difference / — /,„, i.e.. 



exp 



b\I-Ir, 



LM 



(26) 
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FIG. 8. (a) Mean relocation time for different final volt- 
ages, (b) Logarithm of the mean delay time vs current dif- 
ference between final current and the maximum or minimum 
current Im of the initial branch. 

Fig. |^(a), shows the mean relocation time obtained in 
our simulations as a function of Vf. As the final voltage 
approaches that corresponding to Im, the relocation time 
increases. Fig. ^(b) depicts the mean relocation time as a 
function of {I~Im) on a semilogarithmic scale for Vf val- 
ues between 0.737^ and 0.735V^. The solid line denotes 
a linear fit to the data points, that agrees with the expo- 
nential law proposed by Luo et al These figures are 
qualitatively similar to the corresponding ones depicted 
from experimental data in Refs. Q and Quantita- 
tive differences are due to the above mentioned discrep- 
ancies in the tunneling current and the shot noise 
amplitude. Now we focus on the distribution of switching 
times. Typically, delay distributions are either close to 
symmetric Gaussians or they are asymmetric, depending 
on how far Vf is from the limit point of the I-V charac- 
teristics. We have fitted our numerical distributions by 
least squares to either a Gaussian density: 



W{t,T,a) 



exp 



2a2 



(27) 



or to a first passage time (FPT) distribution 



We have observed this dependence in our numerical re- 
sults too. The dimensionless constants b and c are 
b = 64.9866 and c = 1.6717. Im = 1-501 mA is the unit 
of current. In Luo et al's experiments 1^, Im = 136/iA 
(approximately the height of the first maximum of the 
current in the inset of Fig. 1), b = 10.74 (6 times smaller 
than the numerically calculated value) and c = 3.34 (2 
times larger than the numerically calculated value). We 
thus confirm the exponential dependence of the reloca- 
tion time on the current difference and observe a good 
qualitative agreement between numerically and experi- 
mentally obtained values. 



W-(t,2/,/3)dt=^/y^exp (-^\dz. 



where 



Vexp(2/3i) - 1 



(28) 



(29) 



The parameters of these distributions are t (mean relo- 
cation time) and cr (standard deviation) for the Gaussian 
and y and (3 for the FPT distribution. The results of our 
fitting are depicted in Figures || and |l^. 
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IV. CONCLUSIONS 




Time (ns) 

FIG. 9. Time delay distribution for Vf = 0.737 V. Data 
from numerical simulations have been fitted to a FPT distri- 
bution. 



We have studied how the shot noise due to charge 
quantization affects the relocation time of electric field 
domains after a suddent switch of the voltage. We find 
that the mean relocation time depends exponentially on 
the difference between the value of the current at the fi- 
nal voltage and the value of the current at the end of the 
branch corresponding to the initial voltage. The distribu- 
tion function of delay times after a voltage switch changes 
from Gaussian to a FPT distribution as the final voltage 
approaches the limit point of the stationary I-V charac- 
teristics. These results are in qualitative agreement with 
experiments. 
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APPENDIX A: NUMERICAL SCHEME 
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FIG. 10. Time delay distribution for Vf = 0.75 V. Data 
from numerical simulations have been fitted to a Gaussian 
distribution. 

These results agree qualitatively with the experimental 
ones of Rogozia et al's ^ . As in Ref . |^ , our Figures ^ 
and |l^ show that for values of the voltage far away from 
the current jump the time delay distribution changes 
from an asymmetric FPT distribution to a very narrow 
symmetric Gaussian distribution as Vf departs from the 
voltage corresponding to the current jump. These fea- 
tures have a numerical expression in terms of descriptive 
statistics like the mean, the standard deviation or the 
skewness coefficient as shown in the Tables of Appendix 
The numerically calculated largest and smallest delay 
times are also presented. 



This Appendix is devoted to explain some technical de- 
tails of the simulations. The Platen second order scheme 
gives the vector field E"+^ at discrete time t + At as the 
following function of E" at discrete time t Mm : 



2 V \ dt 

N+l 



H(E",^ n At 



N+l 



J (S^ {M\ ) + (Mi ) + 2S^ (E" )) AW^ 



+ J2 (S^'(U;) + S^'(U1)-2S^'(E")) AVF^' 
- J2 (s^(Mi) - S^(Mi)) {{AW^r - At} 



N+l 

+ J2 (S^'(U;)-S^'(U1)) {AW^JAW^'- + K.j} 

r — l,r^j 

Here SJ(-) is the j-th column of S{-), U± E" ± 
S(E")^ \/At, and H and S are evaluated at 



T = E" + H ( E", ^ ) At + ^ S(E")^ AFF^ 

i=i 



N+l 



Mi ^ E" H E 



dt 



At±S^'(E")yAt. 



AW-' are independent gaussian random variables dis- 
tributed with zero mean and variance At, whereas the 
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are independent two point random variables that 



satisfy 



numerical effects (larger in the Heun scheme). The use 
of a numerical scheme that reduces these effects is then 
clearly justified. 



We have used a time step of At — fO^^ (in dimen- 
sionless units) of the same order as the noise amplitude 
a. The values of the random variables V and W have 
been generated through a random number generator im- 
proved by using a seed selector depending on the com- 
puter clock and an algorithm which allows to avoid the 
sequential correlation usual in this sort of generators . 
The Platen scheme is second-order weakly convergent in 
the following sense. Let (7(E) be any sufficiently smooth 
scalar function (with 2(/?-f 1) continuous derivatives pro- 
vided (3 is the order of the scheme). Let us fix the time 
instant at t corresponding to discrete time n. Then 

|(.g(E"))-(.g(E))|<C•(A^)^ 

for any S (0, (5o), where C and i5o are positive con- 
stants. The Platen numerical scheme is certainly more 
complicated and costly than even a stochastic Heun 
(modified Euler) first order scheme. We have had to use 
it to minimize the effects of numerical noise coming from 
floating-point arithmetic (even our high-precision 64-bit 
arithmetic) and that inherent in interpolating our trans- 
port coefficients and contact functions in the boundary 
conditions. In fact, in the absence of the noise, both the 
Heun and the Platen schemes become the well-known 
deterministic Heun (improved Euler) scheme, that is a 
second-order Runge-Kutta method: 



E'^ 



E" 



A< 



[H(E'' 



-H(E" + H(E")Ai)] 



However both schemes differ in their treatment of the 
noise: the stochastic Heun method is weakly first order 
whereas the Platen scheme is second order. The result 
obtained by using the Platen scheme exhibits less disper- 
sion than that reached by the Heun method, as shown 
in Tables | and An appropriate treatment of the 
noise term avoids the presence of artificial numerical ef- 
fects. The effects of the numerical perturbations can be 
illustrated as follows. Let us use the deterministic Heun 
scheme with random initial conditions corresponding to 
disturbances of the stationary field profile at voltage Vb 
and suddenly switch to voltage Vf. The domain reloca- 
tion times have been measured and they give rise to the 
distributions of Figures ^ and [ij. We have compared the 
mean, standard deviation and skewness coefficient of 
these distributions to those corresponding to the use of 
the stochastic Heun and Platen schemes; see Tables |, g 
and III. Notice that the mean relocation times are sim- 



ilar, while the numerical viscosity contributes to scatter 
the results. The shot noise does not change the mean 
values given by the deterministic model, but the disper- 
sion measured by the standard deviation increases due to 
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FIG. 11. Time delay distribution for Vf 
ferent initial conditions. 
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FIG. 12. Time delay distribution for V> = 0.75 for differ- 
ent initial conditions. 



TABLE I. Descriptive statistics of the relocation time dis- 
tributions obtained with the Heun scheme. 



Heun 



Vf = 0.737 



Vf = 0.75 



Lower Limit (ns) 
Upper Limit (ns) 
Mean (ns) 
Standard Deviation 
Skewness coeff. 



77.692 
128.048 

87.863 
6.803 
1.840 



43.775 
45.633 
44.564 
0.299 
0.131 
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TABLE II. Descriptive 


statistics of the relocation time 


distributions obtained with the Platen scheme. 




Platen 


Vf = 0.737 


Vf = 0.75 


Lower Limit (ns) 


77.827 


43.942 


Upper Limit (ns) 


115.025 


44.960 


Mean (ns) 


87.635 


44.541 


Standard Deviation (ns) 


6.339 


0.167 


Skewness coeff. 


1.4237 


-0.2791 


TABLE III. Descriptive 


statistics of the relocation time 


distributions obtained with pertubed initial conditions. 




Vf = 0.737 


Vf = 0.75 


Mean (ns) 


88.916 


44.579 


Standard Deviation (ns) 


1.773 


0.045 


Skewness coeff. 


0.912 


0.255 
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